Identifying the quantum correlations in light-harvesting complexes 

Kamil Bradler and Mark M. Wilde 
School of Computer Science, McGill University, Montreal, Quebec, Canada H3A 2A1 

Sai Vinjanampathy 

Department of Physics and Astronomy, Louisiana State University, Baton Rouge, Louisiana, USA 70803 

Dmitry B. Uskov 
Department of Physics and Engineering Physics, 
Tulane University, New Orleans, Louisiana, USA 70118 
(Dated: July 22, 2010) 

One of the major efforts in the quantum biological program is to subject biological systems to 
standard tests or measures of quantumness. These tests and measures should elucidate if non-trivial 
quantum effects may be present in biological systems. Two such measures of quantum correlations 
are the quantum discord and the relative entropy of entanglement. Ifere, we show that the relative 
entropy of entanglement admits a simple analytic form when dynamics and accessible degrees of 
freedom are restricted to a zero- and single-excitation subspace. We also simulate and calculate the 
amount of quantum discord that is present in the Fenna-Matthews-Olson protein complex during 
the transfer of an excitation from a chlorosome antenna to a reaction center. We find that the 
single-excitation quantum discord and relative entropy of entanglement are equal for all of our 
numerical simulations, but a proof of their general equality for this setting evades us for now. Also, 
some of our simulations demonstrate that the relative entropy of entanglement without the single- 
excitation restriction is much lower than the quantum discord. The first picosecond of dynamics is 
the relevant timescale for the transfer of the excitation, according to some sources in the literature. 
Our simulation results indicate that quantum correlations contribute a significant fraction of the 
total correlation during this first picosecond in many cases, at both cryogenic and physiological 
temperature. 
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I. INTRODUCTION 

Quantum biology aims to understand if, how, and why biological systems exploit quantum-mechanical effects for 
their functionality or for an evolutionary advantage [IHl]. Exemplary biological systems that may exploit quantum 
effects vary from photosynthetic light-harvesting complexes [S], to the avian compass for bird navigation [B], to the 
olfactory system 7^. Ongoing theoretical research indicates that light-harvesting complexes exploit an environment- 
assisted quantum- walk like effect to enhance energy transport [SUllj. the avian compass exploits a radical ion-pair 
mechanism for increased sensitivity of the earth's geomagnetic field p!2Ul4j . and the olfactory system may exploit 
phonon-assisted tunneling for enhanced detection of smell [7] . 

Light-harvesting complexes seem particularly suitable as biological systems to harness quantum-mechanical effects. 
Their lengthscales and energyscales are on the order where we would expect quantum-mechanical laws to apply |15j . 
but what remains less clear is if they can still harness quantum effects such as entanglement even at physiological 
temperature. A recent numerical study addresses this question by showing that light-harvesting complexes could 
demonstrate stronger-than-classical temporal correlations, even at physiological temperature [16] . However, it remains 
an open task to devise and conduct a realistic experimental protocol that demonstrates an irrevocable test of non- 
classical temporal correlations for light-harvesting complexes. 

The aim of the study in Ref . fTF| , as well as that in Refs. [T71 [TB] , is to address one of the important "quantum 
biological questions": Does the biological system exhibit "quantumness" according to a standard test or measure? 
Theoretical machinery from quantum information science |19j . developed specifically for understanding quantum 
computational and communication devices, should be of immense utility in answering this question. Several articles 
have already begun exploiting such tools. Rcf. [16 exploits the Leggett-Garg test of non-classicality [5D] to suggest 
that light-harvesting complexes might exhibit stronger-than-classical temporal correlations, while Refs. [T71IIH] utilize 
standard entanglement-based measures of quantum correlations [2TH23] to suggest that they might exhibit stronger- 
than-classical spatial correlations. Such standard measures of quantum behavior are more convincing than, say, a 
claim that wavelike motion in population elements of a density matrix is a signature of quantumness |24) . 

The studies in Refs. [T71 [TH] both suggest that it might be possible to observe spatial quantum correlations in 
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the Fenna-Matthews-Olson (FMO) light-harvesting protein complex, but the quantum correlation measures exploited 
there, such as concurrence |25j . the measure based on the global relative entropy of entanglement [23], and logarithmic 
negativity , might not capture all of the spatial quantum correlations that are present in a given quantum system. 
For example, consider a bipartite quantum system in the state: 

^|O)(O|^0|+>(+|^ + ^|->(-|-'0|l)(l|^, (1) 

where |±) = (|0) ± |1)) /\/2 and Alice possesses the system A and Bob possesses the system B. The state in ([l]) is a 
separable state |27| , meaning that Alice and Bob can prepare it by means of local quantum operations and classical 
communication. Yet, it is fundamentally non-classical because the states |0) and |— ) on Alice's local system or |1) 
and |-|-) on Bob's local system are indistinguishable from one another (they are non-orthogonal). It could be possible 
for two bacteriochlorophyll sites of a light-harvesting complex to admit a state of the above class or a mixture of such 
states, but the quantum correlation measures mentioned above all vanish for such a separable state, despite its status 
as a fundamentally non-classical state. Thus, these measures might not capture the full "quantumness" that might 
be present in a light-harvesting complex. 

A different measure of quantum correlations, known as the quantum discord 28, 29 , captures all correlations in a 
quantum state that are non-classical. The discord is a measure of the total correlations present in a shared quantum 
state, reduced by the classical correlations obtainable when one party performs a local measurement. For example, the 
state in ([T]) registers a non- vanishing quantum discord and therefore possesses non-classical correlations according to 
this measure, even though it does not violate a Bell inequality 30J due to it being a separable state. Other examples 
demonstrate that the quantum discord in a quantum system can remain positive even if the entanglement vanishes 
after a finite time [3T]. States that register non- vanishing discord can lead to exponential speedups [32 j^S] in a 
quantum computational model known as the "one-clean qubit model" |34j , where the only requirement for a speedup 
is access to a single qubit in a pure state. Although it is unlikely that a light-harvesting complex could be performing 
exotic quantum computational speedups of the aforementioned nature or in the standard way |35J , one cannot rule out 
the possibility that a light-harvesting complex exhibiting non-vanishing discord could be exploiting "quantumness" 
of this form for enhanced energy transport. 

In this paper, we systematically study the presence of quantum correlations in the FMO light-harvesting complex, 
both at cryogenic and physiological temperature. Our first contribution is a simple formula for computing the relative 
entropy of entanglement [SIT when a light-harvesting complex evolves according to the open quantum system model 
well studied in the quantum biological literature llj. This formula assumes that dynamics are restricted to a zero- 
and single-excitation subspace, and it reduces the a priori computationally intensive optimization task for computing 
the single-excitation relative entropy of entanglement to a simple calculation with quantum entropies. These results 
for the single-excitation relative entropy of entanglement generalize those of Sarovar et al. in Ref. [T7| for the global 
relative entropy of entanglement. We then calculate the quantum discord for several phenomenologically motivated 
"bipartite cuts" of the FMO protein complex and find that it is equivalent to the single-excitation relative entropy of 
entanglement in all cases presented here. These results demonstrate that quantum correlations contribute a significant 
fraction of the total correlation during the highly relevant first picosecond of dynamics in many cases (some sources 
in the literature [101 111] indicate that the average time it takes for an excitation to trap to the reaction center is one 
picosecond after it arrives from the chlorosome antenna). In other cases, these quantities are equivalent and contribute 
a non-negligible fraction of the total correlation during the first picosecond. This contribution indicates that non- 
classical spatial correlations may be playing a role in the efficient transfer of an excitation from the chlorosome antenna 
to the reaction center. Our final contribution is to study the relative entropy of entanglement without its optimization 
restricted to the single-excitation subspace, and we find that it can be significantly less than the restricted relative 
entropy of entanglement while only using just a small fraction of doubly-excited states in the optimization. 

We structure this paper as follows. We first briefiy review both the dynamical model of the FMO protein complex, 
the definition of the quantum discord, the definition of the relative entropy of entanglement, and we then discuss 
how to compute these quantities in the FMO protein complex model. Our analytic result states that the relative 
entropy of entanglement, when restricted to the zero- and single-excitation subspace, is equal to a simple formula 
that is a difference of entropies. Appendix [A] provides a proof of this theorem. Section [TV] p resents the results of our 
numerical simulations under various initial configurations and temperatures, and Section^ discusses the unrestricted 
optimization of the relative entropy of entanglement. We then conclude with observations and open questions for 
future research. 



3 



II. REVIEW 
A. FMO Complex Dynamics 

The FMO protein complex is the crucial, light-harvesting component of the green sulfur bacteria prosthecochloris 
aestuarii, that develop in dimly-lit, anoxic environments such as stratified lakes or sulfur springs |37) . It is a trimer, 
consisting of three identical subunits. We study one unit of the trimer, consisting of seven bacteriochlorophyll sites 
that act as a "molecular wire," transferring energetic excitations from a photon-receiving antenna to a reaction center. 
A photon impinges on the antenna, producing an electronic excitation, dubbed an exciton, that then proceeds to the 
unit with seven sites. While traversing the seven sites, the exciton can either recombine, corresponding to an energetic 
loss, or it can trap to the reaction center for energy storage. 

We characterize the quantum state of the exciton as a density operator in the site basis: 

m.nG{G,l.,...,7.,S} 

where the state \m) indicates that the exciton is present at site m. The sites can be any of the seven sites in the 
protein complex m G {1, . . . , 7}, a ground state \G) that represents the loss or recombination of the exciton, or a sink 
state 15) that implies that the exciton has trapped to the reaction center. We adopt the following "qubit" convention 
in this work (as adopted in previous works [T71 [SHI El] ) , by assigning site states to tensor-product states: 

\G) = 15)2 19)3 15)4 13)5 15)5 Iff)? \9)s ' 
|1) = |e)i 15)2 1^)3 1^)4 lff)5 lff)6 Iff)? Iff)s ' 



|7) = |g)i 15)2 lff)3 lff)4 lff)5 lff)6 |e)? Iff)s ' 

\S) = \g)i 15)2 lff)3 lff)4 lff)5 lff)6 Iff)? \^)s ' 

where g indicates the absence of an excitation and e indicates the presence of an excitation at a particular site. The 
excitation number is a conserved quantity in the absence of light-matter interaction events |40j , and this observation 
restricts the protein dynamics to a zero- and single-excitation subspace. Thus, the above nine states are the only 
states that we consider for the exciton while it traverses the seven sites. 

Observe that "tracing over any site" in this qubit representation has the effect of placing the population term for 
that site into the population of the ground state. For example, suppose that we trace over all sites except for the first 
one. The resulting density matrix has the form: 

{pGG + P22 + . . . + P?7 + Pss) \G) (G| + pn |1> (1| + piG |1) (G| + pGi \G) (1| . 

Such manipulations are important for computing correlation measures for any bipartite cut of the sites in the FMO 
complex. 

Evolution of the density matrix occurs according to a combination of both coherent and incoherent dynamics. A 
tight-binding Hamiltonian of the following form governs coherent evolution across the seven site states . . . , |7): 

H = ^E^ |m) (to| -I- Vnm {\n) (m| + \m) {n\) , 

m n<m 

where is the relative energy at site m and Vnm represents a coupling between sites n and m (see Refs. [TOl [HI [H] 
for the Hamiltonian that governs evolution). The third site has the lowest energy of all seven sites and is closest 
to the reaction center, and it is for these reasons that researchers suggest that the objective of the molecular wire 
is to transmit the excitation as quickly as possible to the third site [SHTTj. The off-diagonal terms lead to coherent 
couplings between sites — the strongest couplings are between the first and second sites, the fourth and fifth sites, and 
the fifth and sixth sites. Recent work determines the energy landscape of the FMO complex that reveals coherent 
pathways the excitation may traverse to get to the third site [231 [S] ■ 

A set of three Lindblad superoperators governs incoherent evolution of the FMO complex, and each superoperator 
has a phenomenological justification for its presence in the evolution [SV HOl ITS] . The first incoherent mechanism iZdiss 
is due to recombination or loss of the exciton, acting like an amplitude damping mechanism. The second incoherent 
mechanism >Csink is due to trapping of the excitation from the third site to the reaction center. The final decoherence 
mechanism £doph is a local dephasing at each site, due to unavoidable interaction with the surrounding protein 
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environment. Ref. [3] shows how to relate the dephasing rate to the temperature of the system, under a particular 
decoherence model (we employ this relation in our numerical study in this paper). Refs. [8^1101 [16] describe in detail 
the decoherence mechanisms in the FMO complex. Despite the Markovian nature of the evolution, it describes the 
dynamics reasonably accurately [S] (though see the recent work on non-Markovian evolution in the FMO complex 

mm)- 

Putting everything together, the following master equation gives the full evolution of the density matrix p: 

p= -i [H, p] + £diss (p) + >Csink (p) + >Cdoph (p) ■ (2) 

Note that the above evolution induces coherences only between sites in the FMO complex, so that any off-diagonal 
element of the density matrix of the form pi^c or /9i,s, where i € {1, . . . , 7}, vanishes. Various sources in the literature 
[9] and recent experimental evidence [43] suggest that the initial state of the FMO complex is a pure excitation at 
site one, a pure excitation at site six, or a mixture of these two states because sites one and six are closest to the 
antenna. We incorporate each of these cases in our simulations. 

The dynamical model in ^ implies that the FMO complex exhibits both classical and quantum random walk 
behavior, but does not exhibit exotic "quantum stochastic walk" behavior according to Table I of Ref. [44 . The 
classical and quantum random walk behavior then suggests that this evolution should establish both classical and 
quantum correlation between sites in the complex. It is these correlations that we measure in this paper. 



B. Quantum Discord 

We first recall various informational measures of a quantum state before briefly reviewing the motivation for the 
quantum discord as a measure of quantum correlations. Suppose that two parties Alice and Bob share a quantum 
state p^^ . The von Neumann entropy of this state is as follows: 

H{AB)^^-tr{p^^\ogp^^}, 

where the logarithm is base two. The entropy H {AB)^ measures the uncertainty about the total quantum state p^^ 
in units of bits, whenever the logarithm is base two |45j . Similarly, we can compute the marginal entropies H (A) and 
H (B) — these entropies are with respect to the respective reduced states p"^ and p^ , obtained by a respective partial 
trace over Bob or Alice's share of the state p'^^ . 

The quantum mutual information / ( A; B) is a measure of the total correlations, both classical and quantum, shared 
between two parties, where 

I{A-B)^^H{A)^ + H(B)^-H{AB)^. (3) 

The quantum mutual information admits a natural operational interpretation as the amount of noise necessary to 
destroy all correlations present in a bipartite state |46) . Its other operational interpretations are as the maximum 
amount of secret information that one party can send to another if they use a shared state as the basis for a one-time 
pad cryptosystem |47j, the entanglement-assisted classical capacity of a quantum channel |48H50j . and more recently 
the amount of private quantum information that a sender can transmit to a receiver while being eavesdropped on by 
a uniformly accelerating third party [51]. All of these operational interpretations justify the notion of the quantum 
mutual information measuring correlations between two parties in units of bits. 

Suppose now that Alice would like to extract classical information from the quantum state p"^^ . She performs a 
von Neumann measurement {|a;) on her share of the state, where the states {|a;)} form an orthonormal basis. If 
Alice obtains classical result x from the measurement, the resulting conditional quantum state of the whole system is 

\x) (xl^^pf, 

where 

^ -^Trx U\x) {xf ® /^)p^^(|x) {xf (4) 
p [x] I J 

p(x)^Tr{|x)(x|^p^}, (5) 

and we now label Alice's system with X because it is classical. If we then take the expectation over all of Alice's 
outcomes, the description of the system is a classical-quantum state: 

a^^EE^p(x)|a;)(xr®pf. (6) 

X 
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We might think that a natural measure of the classical correlations present in a bipartite state is the amount of 
correlations in the resulting classical-quantum state: 

The above expression features a maximization over all of the von Neumann measurements A"* that Alice could perform. 
These measurements result in a state of the form a in (j6|, and the quantum mutual information I {X;B)^ is with 
respect to such a state. In fact, Henderson and Vedral proposed such a measure [3S], and Devetak and Winter later 
justified this measure by providing an operational interpretation of it as the maximum amount of common randomness 
(perfect classical correlations) that two parties can extract from a bipartite quantum state |52) . 

The two correlation measures in ([s]) and ([7|), the first a measure of total correlations and the second a measure 
of classical correlations, suggests that a measure of the quantum correlations should be the difference of these two 
quantities. The quantum discord D{A;B)p is such a measure (SUES], defined as the difference of total and classical 
correlations: 



D{A;B)^^I{A;B)^-I,{A;B)^. (8) 

The discord is always non-negative, by the quantum data processing inequality of quantum information theory |54) . 
and it is generally not symmetric with respect to the parties A and B: 

3p:D{A-B)^^D{B:A)^. 

It captures all of the quantum correlations, including entanglement, but also the quantumness in separable states of the 
form in ([T]) . Zurek has suggested a physical interpretation of the discord as the difference in efficiency between quantum 
and classical Maxwell's demons [SS], but it still lacks a clear operational interpretation in the sense of Refs. 
Nevertheless, we still employ the discord as a measure of the quantum correlations in the FMO complex. 
We can rewrite the quantum discord as the following expression: 

D {A- B)^ = I {B)A)^ + min (a:) H (pf ) , (9) 

X 

by performing straightforward manipulations of the entropies in ([s]). In the above, / {A)B) is the coherent informa- 
tion [SI], equal to the negative of a conditional entropy: 

/ {B)A)^ = -H {B\A)^ = H {A)^ - H [AB)^ , 
and the probabilities p {x) and conditional density operators are as they appear respectively in Q and ([5]). 



C. Relative Entropy of Entanglement 

The relative entropy of entanglement is an entanglement measure from quantum information theory f36', and we 
briefly review its definition. The relative entropy D [p^^\\uj^^^ of two bipartite states p"^^ and uj"^^ is as follows |54j : 

D {p^^Wu^^) ^ Tr {/^ logp^^} - Tr {/^ logc.^^} . (10) 

This measure, in some sense, quantifies the "distance" between two bipartite states, but it is not a distance measure 
in the strict mathematical sense because it fails to be symmetric. Though, this intuition is useful for constructing 
an entanglement measure. We might naturally expect a good measure of entanglement to be the distance of a given 
bipartite state to the closest separable state = '^iP{i) erf (8 af . The relative entropy of entanglement R (p^^) 
is such a measure, defined as the minimization of the relative entropy over all separable states: 

i?(0^mini?(/^||a-4^), 

where S is the class of separable states. This measure satisfies the properties of an entanglement measure and appears 
extensively in the quantum information theory literature [36J. We mention that Ref. [T7^ employed this entanglement 
measure for determining quantum correlations in the FMO complex, but they considered the global relative entropy 
of entanglement rather than a bipartite version of it as we consider in this work. 
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III. FORMULA FOR THE SINGLE-EXCITATION RELATIVE ENTROPY OF ENTANGLEMENT 

The restriction of dynamics to the zero- and single-excitation subspace allows for significant simplifications to 
the theory of energetic transfer in the FMO complex. We can also apply this restriction to the relative entropy of 
entanglement, by restricting the optimization over separable states to the zero- and single-excitation subspace. Let 
Re (p^^) denote the relative entropy of entanglement with this restriction applied to its optimization. Theorem [l] 
below states that Re {p"^^) is equal to a simple formula that is a difference of entropies. Thus, the theorem significantly 
simplifies the computation of this quantity. 

Theorem 1 Consider a density operator p^^ restricted to the zero- and single- excitation subspace. Let 

A ^ a \G) (G| +p^® \G){Gf + \G){G\^ ® , 

where a is the population of the ground state, p^ is the projection of Alice 's part of p^^ into the single- excitation 
subspace, and pf is defined in a similar way. Suppose that the dynamics for this density operator never induce any 
coherences between the zero- and single- excitation subspaces (as in the dynamics in Then the single- excitation 

relative entropy of entanglement Re (p) is equal to the difference of the entropy of A (p^^) and p"^^ : 

Re 

The full proof appears in Appendix [Xj It exploits standard properties of the von Neumann entropy and a pertur- 
bative argument. 

IV. SIMULATION RESULTS 

We conducted several simulations at both cryogenic temperature (77°K) and physiological temperature (300°K) and 
for the initial state being a pure state at site one, six, and the mixture of the previous two. These simulations calculate 
both the quantum mutual information, quantum discord, and single-excitation relative entropy of entanglement with 
respect to several "bipartite cuts" of the sites in the FMO complex. Throughout this section, we refer to quantum 
discord, but all our simulations indicated that the single-excitation relative entropy of entanglement is equal to the 
quantum discord. These results provide strong evidence that the quantum discord is equal to the formula in Theorem[T] 
for these cases, but a general proof of this conjecture eludes us for now. Our different simulation cases are as follows: 

1. We first considered the cut where system A consists of site three and system B consists of sites one and six. We 
picked this cut because the initial state of the complex is at site one, six, or the mixture, and the objective of 
the FMO "molecular wire" is to transfer the excitation from these initial sites to site three. If spatial quantum 
correlations play a role in the transfer of the excitation, one would expect a state with this bipartite cut to 
register a non-negligible amount of quantum discord. 

2. The next bipartite cut that we considered is with the A system consisting of sites one and two and the B system 
consisting of site three. We picked such a cut because recent analysis of the FMO Hamiltonian 4f suggests 
that a superposition state of sites one and two gives an efficient energetic pathway for the excitation to transfer 
to site three. The energy of the superposed state is closer to the energy of site three than it is to the energy of 
either site one or site two. We would again expect such a cut would register a non-negligible amount of quantum 
discord if quantum correlations play a role in the transfer of the excitation. 

3. Finally, after conducting the above simulations, we conducted simulations with respect to the cut where system 
A consists of site three and system B consists of all other sites. The quantum mutual information for this case 
should be larger than for the case of the other cuts, by the quantum data processing inequality. 

Figure [T] plots the results of the first simulation for cryogenic temperature and for each initial state. The figure 
indicates that the quantum discord contributes a significant fraction of the total correlation for short timescales (less 
than one picosecond). This contribution of the quantum discord to total correlation is significant, considering that 
the average transfer time of the excitation to site three occurs around the order of one picosecond in our model. For 
longer timescales (greater than one picosecond) , the total correlation between A (sites one and six) and B (site three) 
increases to its maximum at around 2-4 ps, and the quantum discord no longer contributes any significant amount to 
the total correlation. Even though the total correlation rises so much higher for longer timescales than it is for shorter 
timescales, the increase and peak are not relevant for excitation transfer, i.e., "they arrive too late" given that this 
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FIG. 1: (Color online) Mutual information I {A; B) and quantum discord D {A; B) with system A as site three and system B 
as sites one and six. The figure plots these quantities at cryogenic temperature (77° K) as a function of time when the initial 
state is (a) a pure state at the first site, (b) a pure state at the sixth state, (c) an equal mixture of the two previous states. In 
each of the above plots, the quantum discord is a significant fraction of the total correlations during the first picosecond. 
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FIG. 2: (Color online) Mutual information I {A; B) and quantum discord D [A; B) with system A as site three and system B as 
sites one and six. The figure plots these quantities at physiological temperature (300° K) as a function of time when the initial 
state is (a) a pure state at the first site, (b) a pure state at the sixth state, (c) an equal mixture of the two previous states. In 
each of the above figures, the quantum discord is a non-negligible fraction of the total correlation during the first picosecond. 




FIG. 3: (Color online) Mutual information / [A; B) and quantum discord D {A; B) with system A as site three and system B 
as sites one and two. The figure plots these quantities at cryogenic temperature (77° K) as a function of time when the initial 
state is (a) a pure state at the first site, (b) a pure state at the sixth state, (c) an equal mixture of the two previous states. 
In each of the above cases, the quantum discord is a significant fraction of the total correlation during the first picosecond. In 
fact, for the second case, the quantum discord contributes nearly all of the total correlation during the first picosecond. 
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FIG. 4: (Color online) Mutual information I (A; B) and quantum discord D {A; B) with system A as site three and system B as 
sites one and two. The figure plots these quantities at physiological temperature (300°K) as a function of time when the initial 
state is (a) a pure state at the first site, (b) a pure state at the sixth state, (c) an equal mixture of the two previous states. In 
each of the above cases, the quantum discord contributes a fraction of the total correlation during the first picosecond. In the 
second case, it again contributes nearly all of the total correlation during the first picosecond. 



transfer occurs on the order of 1 ps. The total correlation then washes away as time increases beyond 10 ps, and it 
should nearly vanish for times around 1 ns because this time is the average recombination time of the exciton, and 
no correlations should persist after a recombination. 

Figure [2] plots the results of the first simulation for physiological temperature. The difference between cryogenic 
and physiological temperature is a qualitative shrinking by a factor of two, but the fraction of quantum discord that 
contributes to total correlation for short timescales is lower than it is for cryogenic temperatures. The local dephasing 
at each site acts to destroy both total and quantum correlation, but it appears to have a more harmful effect on 
quantum correlation. Despite the low amount of quantum discord registered at physiological temperature, it could 
still be that a light-harvesting complex is harnessing this small amount of quantumness (that is a significant fraction 
of total correlation) on this short timescale to enhance excitation transfer. 

Figures [3] and |4] plot the results of the second simulation for both respective temperatures and for each initial 
state. These results are qualitatively similar to the previous bipartite cut, but the most striking difference is that the 
quantum discord contributes all of the correlation (it is equal to the quantum mutual information) for short timescales 
when the initial state is at site six for both cryogenic and physiological temperatures (see Figures |3][b) and[4]Jb)). 

Figures [5] and |6] plot the results of our final simulation, where the bipartite cut has system B as site three and 
system A as all other sites. The results for quantum discord are again qualitatively similar to previous results, but 
the quantum mutual information is significantly higher than for the previous cases (this is expected because the 
quantum data processing inequality states that correlations can only decrease when discarding subsystems). The 
quantum discord again contributes a significant fraction of the total correlation for short timescales and contributes 
only a small fraction for later timescales. Again, physiological temperature decohcrcnce mitigates the presence of 
quantumness. 



COMPARISON OF THE QUANTUM DISCORD WITH THE STANDARD RELATIVE ENTROPY 

OF ENTANGLEMENT 



In this section, we focus on computing the relative entropy of entanglement in ( 10 1 without the assumption that the 
set of separable states involved in the minimization of D[p'^^\\<7'^^) is limited to the single-excitation subspace. 
This implies that may include multiply-excited states (for example, states for which both systems A and B carry 
an exciton). 

A question may arise as to why such multiply-excited states are relevant for the optimization of the relative entropy 
of entanglement. After all, the Hamiltonian of the system commutes with the operator corresponding to the total 
number of excitations in the system, and the dissipation operators in ^ only decrease the number of excitations 
by dumping them to the sink or to the reaction center. There are two different ways to answer this question. The 
first one is simple yet somewhat formal — Vedral established the relative entropy of entanglement in the context of 
quantum information theory as a measure of entanglement for a general set of states with no subspace limitations 
imposed. Therefore, even though the form of specific Hamiltonian and dissipation terms in ([2| leads to the absence of 
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FIG. 5: (Color online) Mutual information I (A; S) and quantum discord D (A; B) with system A as site three and system B 
as all other sites (1-2, 4-7). The figure plots these quantities at cryogenic temperature (77°K) as a function of time when the 
initial state is (a) a pure state at the first site, (b) a pure state at the sixth state, (c) an equal mixture of the two previous 
states. The quantum discord contributes a significant fraction of the total correlation during the first picosecond. 
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FIG. 6: (Color online) Mutual information / [A; B) and quantum discord D [A; B) with system A as site three and system B 
as all other sites (1-2, 4-7). The figure plots these quantities at physiological temperature (300°K) as a function of time when 
the initial state is (a) a pure state at the first site, (b) a pure state at the sixth state, (c) an equal mixture of the two previous 
states. The quantum discord does not contribute a large fraction of the total correlation here. 



multiply-excited states in the density operator, there is no particular reason to consider that multiply-excited states 
and the mode of preparation of a state are irrelevant in this setting. From this perspective, entanglement is determined 
solely by the state rather than by the history and mode of state evolution. The relative entropy of entanglement is a 
measure of the "distance" of a state to the closest separable state and excluding multiply-excited states from the set 
of separable states introduces a distortion to the original concept. Actually, the measure of such a distortion may be 
rather large, as we found out performing optimization of D{p^^\\a^^) in the full space (see Figure [rj^ a)). 

Our second answer is more physically motivated. In general, while analyzing the exciton dynamics in a light 
harvesting system, one assumes that the intensity of light, which determines the photon flux, is quite small. Then it 
follows that the probability of simultaneous creation of two excitations is negligible. By performing the optimization 
of the D{p^^\\(j^^) in the full space, which includes multiply-excited states that can have the form \E)^ , and 
computing the amount of multiply-excited state population in the closest separable state tr^^, we can set up a limit 
for the light intensity when the single-excitation assumption will become meaningless even in the context of excitation- 
preserving dynamics. In accordance with Caratheodory's theorem, we generated the generic separable state as a 

64 

sum of 2^ = 64 arbitrary pure separable states ^ c„ IV'n) (V'nl- Here, we limited our computation to a simple example 

n=l 

when site A — \,Q and site 5 = 3, such that 



\-^n) = [ao |5)6 + "1 |5)6 + "2 \g)i |e)6 + "3 |e)i |e)g) ® (/3i \g)^ + P2 \e)^) 
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FIG. 7: (a) This figure simulates the dynamics of the FMO complex at 77°K and calculates both the single-excitation relative 
entropy of entanglement and the full relative entropy of entanglement for A = 1,6 and B = 3. (b) The ratio of the number 
of doubly-excited states to the number of singly-excited states in the optimal separable state for the relative entropy of 
entanglement is rather small. 



The optimization is performed by the gradient method in the space of coefficients c„, a^, and /3j. Figure [7|[a) shows 
results of the computation for T — 77° K. We note that relative entropy of entanglement for the single-excitation- 
subspace case has a similar shape to the full-space relative entropy of entanglement, but it is approximately five 
times smaller. Figure [Tf^b) demonstrates that the population of doubly-excited states in the optimal separable state is 
negligibly small — it is three to five orders of magnitude smaller than the population of singly-excited states. We find 
it remarkable that such a negligible admixture of doubly-excited states changes the relative entropy of entanglement 
by more than a factor of five. 

VI. CONCLUSIONS 

We presented results quantifying the quantum correlations present in a biological system at cryogenic and physi- 
ological temperature. Theorem [T] proves that the single-excitation relative entropy of entanglement admits a simple 
form. We then simulated the dynamics of the FMO complex and calculated the quantum discord and the quantum 
mutual information for various phenomenologically motivated bipartite cuts of the sites in the FMO protein complex. 
It is surprising that the quantum discord and the single-excitation relative entropy of entanglement are equivalent 
for these simulations, but we have not been able to find a general proof that the quantum discord of this system is 
equivalent to the formula in Theorem [T] The results of our simulations indicate that quantum correlations contribute 
a significant fraction of the total correlation during the first picosecond of dynamics in many cases. Our last contribu- 
tion was to study the relative entropy of entanglement with an unrestricted optimization, and we found that a small 
fraction of doubly-excited states contribute significantly to reducing the relative entropy of entanglement. 

Open questions remain for this line of research. Our last contribution above suggests an intriguing open question 
of whether optimizing the relative entropy of entanglement over a single-excitation subspace, as done here and in 
Ref. [IT, is operationally justified. We should also compute the correlation of the quantum discord with the measure 
of transfer time in Ref. j9j of the excitation from the antenna to the reaction center. This correlation might indicate 
how relevant "quantumness" is for the transfer of the excitation. Finally, it might be interesting to determine how 
efficiency or transfer time is affected by differing amounts of quantum discord, if there can be efficient energy transfer 
without quantum discord, and how discord compares with entanglement in the general case. 
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Appendix A 



We first establish some notation before proving Theorem [T] Let our zero- and single-excitation space be spanned 
by the following states: 



\G) 



AB 



,AB 



l.9>i 1.9)2- 
|e)i 15)2- 



\9)n ' 



The number of states is thus n + 1. We introduce a bipartite splitting of the above states so that the corresponding 



subsystems A and B are spanned by { | G) 
convention as above. Note that n 



ID 



Knand{|G)^,|l) 



\nb) }, respectively, using the same 



where k = 1 , . . . , ria and j — 1 , . . 
individual subsystems of A and B: 



+ 1 — n, and the original basis written in terms of the split bases reads 
{|G)^|G)^, |fc)^|G)^, |G)-^|j)^}, 
rif). Let and 11^ denote the projectors onto the excited subspaces of the 



n 



k=l 



n 



rib 



|j>Of 



The projector onto the full bipartite single-excitation subspace is as follows: 

n^^ = |G) (G|^ ® nf + ® |G) (G|^. 



(Al) 



(A2) 



Let = |G)(G|"^^. Then 4- is a projector onto the zero- and single-excitation subspace. We have that 
(G|'^^p'^^|fc)"^^ = because the dynamics in ^ do not induce any correlations between the ground state and the 
single-excitation subspace of the density matrix p^^ . The following chain of inequalities then gives another way to 
write the density matrix p^^ that is more useful to us: 



(n 



= aU, 
= aU, 

= an, 



AB 

a 

AB 

a 

AB 

a 



ttAB ABt-iAB 



(|G)(G|^®nf)p^^(|G)(G| 



n^®|G)(G|^)/^(|G)(G| 
ab 



n 



nf + n^®|G)(G|^) 
(|G)(G|^®nf)p^^ (nf ®|G)(G|^) 



+ (n^ ® |G)(G|^) p^^ (|G)(Gr ® nf ) + (n^ ® |g)(g|^) p^^ (nf ® |g)(g| 



an^ + pf ® |G)(G|^ + |G)(GK 



-AB 



(A3) 



The first equality follows because p^^ lives in the zero- and single-excitation subspace. The second equality follows 
by expanding. The third equality follows because the density operator p^^ does not have any correlations between 



the ground state |G) and any excited state. The fourth equality follows from the relation in (A2). The fifth equality 
follows by expanding, and the last follows from the following definitions: 

jA ^ /r^\B\ AB (tjA ^ //-<|S^ 



p--^(nf ®|G)(G|^)p^^ (nf ®|G)(G|^ 



Pe 
AB 



n^ 



(|G)(Gr0nf)p^^(|G)(Gr 

(nf ®|G)(G|^)p^^(|G)(Gr®nf) 
+ (|G)(G|^0nf) p^^ (nf 0|G)(Gr) 



(A4) 
(A5) 
(A6) 
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The matrix pf is the projection of p"^^ into Ahce's single-excitation subspace, is similariy defined, and r^^ is the 
block-off-diagonal part of p^^ living in the single-excitation subspace. 

1. Proof that {p^") = H {A (/^)) - H {p-"'') 

Proof. Our first candidate to maximize the relative entropy of entanglement is a state of the following form 

ao = l3Ug + af® \G){G\^ + \G){G\'^ ® af , (A7) 

because it is in the most general form of a single-excitation separable state with no correlation between the ground 
state and the other states. The operators cr^ and af are some positive operators that live in the single-excitation 
subspaces with respective projectors 11^ and Ilf. So for now, we restrict the minimization in the relative entropy of 
entanglement to states of the above form, leading to the following minimization: 

i?e (p-^^) ^ minTr {p^^ logp^^} - Tr {p^^ logao} . (A8) 
We first determine a spectral decomposition of and as follows: 

Ua fib 

at^Y.ak\4>k){c^k\^. (A9) 

k=l k=l 

Let OJ^ = log (c^) and uj^ = log (erf). The following equalities then hold 

r^^ log ao =[{nf®\G) (G^) p^^ (|G) (G^ ® nf ) + (|G) (G^ ® nf ) p^^ (nf <E> |G) (G|^)] 
log {aUg + (7^ ® |G)(G|^ + |G)(G|-^ ® af ) 
= (nf |G)(G|^) p-^^ (|G)(G|^ ® c^^) + (|G)(G|^ ® nf ) p^^ (w^ \G){Gf) 

The first equality follows by definition, and the second follows by considering the overlapping support of the different 
subspaces. We then find that 

Tr {r^^ log ao} = Tr { (nf « |G)(G|^) p^^ (|G)(G|^ <E> c.^) + (|G)(G|-4 <E> nf ) p^^ {oj^ <g> |G)(G|^) } 

= Tr {(|G)(G|^ c.^) (nf |G)(G|^) p-^^} + Tr {(u;^ ^ |G)(G|^) (|G)(G|^ nf ) p^^} 
= 0. 

The second equality follows from cyclicity of the trace, and the last equality follows because the zero- and single- 
excitation subspaces are orthogonal. Therefore, 

Tr{p-4^1ogao} =Tr{(p^^ -T^^) logao}. (AlO) 

Observe that p"^^ — t^^ is a state of the form ao with (3 = a, = p^, and af = pf . Consider that 

Tr {(p-^^ - r^^) log(p-4^ - r^^)} > Tr {{p^^ - r^^) logao} , 



because the relative entropy D {{p"^^ — '''"^■^)||ao) > [53]. We conclude that Re {p^^) hi (A8) attains its global 
minimum at p — r"*^. Thus, 

Re {p^^)^H{A{p^^))^H{p^^). 



The expression in ( A7 1 is not the most general single-excitation separable state because it does not include coherence 
elements between the ground state and the single-excitation subspace. That is, matrix elements of the following form 
could be nonzero: 

|G)^|G)^(0fe|^(G|^ + H.C. (All) 
|G)^|G)^(G|-^(0fc|^+H.C. (A12) 
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We now show that it is not necessary to consider states o f th is more general form by appeahng to a perturbation 



theory argument [56]. Consider a sohition uj = p — t^^ of iA8h. Let us take a small shift of this solution: 

Lj^uj + e{ \Gf \Gf {cj^kl^ {Gf + H.C.). (A13) 

Consider that 

L.|0,)^|G)^=afe|^fe)^|G)^, (AM) 

and 

o.|G)^|G)^ = a|0fc)^|G)^. (A15) 

Since u is Hermitian, we can apply perturbation theory to compute the new states and energies. Brillouin-Wigner 
perturbation theory states that if such a Schrodinger-like equation above is perturbed via 

Lu ^ Lu + eduj. (A16) 

where duj = \G)^\G)^ {(j^k]"^ {G\^ + H.C., then the new states are given by 

\^k)^\G)'' ^\4>k)^\G)^ +e 1-^™) ^-(0™|dc.|0fc)^|G)^+O(62). 

Here Em refers to the eigenvalue of the state \4>m)- Substituting for dw, we obtain 

Similarly, applying BWPT to the second Schrodinger like equation above yields 

|G)^|G)^^|G)-^|G)^ + e l'/'m>^T^(0™Mc.|G)^|G)^ + O(e2). 



a — En 

|0„>#|G>-*|G>« 



Again, substituting for duj gives 



|G)-4|G)^ ^ |G)^|G)^ + ^^|^fe)-4G)^ + 0{e^). 



Within BWPT, the eigenvalues are calculated as 

CLk — E, 



a,^a, + e(0,|^(G|^dc.|0fe)^|G>^+ V — |(0„,|du;|(/.fc)^|G)^p 

^ — ^ n 1. — 



|0„)#|fc>-*|G) 

and similarly 



2 



a^a + e{G\^{Gfdij\G)^\G)''+ V — |(0,„|du;|G)-^|G)^'2 

^ — ' a — Em 

|0„)#|G>-*|G)-B 

Using the above, we find that the eigenvectors transform up to first order as follows: 

\Gf ^ 10.)^ \Gf + eMlM! + oie% (A17) 

\Gf \Gf ^ |G)^ \Gf + e^^Li^ + oie'). (A18) 

a — ap^ 

Correspondingly, the eigenvalues change up to second order as follows: 

afc->afe + 0(e2), (A19) 
a^a + 0{t^), (A20) 
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and so 



e , a 
w — )• w H log I — 



\Gf\Gf{cl>kf{Gf+^.C. +0{e') 



Using 



Tr \Gf{<Pk\^ {Gf + H.C.)] = 

we conclude that oj = p — t^^ is the state where the relative entropy attains its global minimum. 



(A21) 
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